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ABSTRACT 

In this paper we present a study of the cosmic flow from the point of view of how 
clusterings at different dynamical regimes in an expanding universe evolve according 
to a 'coarse-grained' partitioning of their ranked energy distribution. By analysing a A- 
CDM cosmological simulation from the Virgo Project, we find that cosmic flows evolve 
in an orderly sense, when tracked from their coarse-grained energy cells, even when 
nonlinearities are already developed. We show that it is possible to characterize scaling 
laws for the Pairwise Velocity Distribution in terms of the energy cells, generally valid 
at the linear and nonlinear clustering regimes. 
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1 INTRODUCTION 

Cosmological N-body simulations play a crucial role in the 
study of the formation of cosmic structure. These struc- 
tures are believed to have grown from weak density fluc- 
tuations present in the otherwise homogeneous and rapidly 
expanding early universe. These fluctuations were amplified 
by gravity, eventually tu rning into the rich s tructure that 
we see around us today iSprineel et al.ll2005l) . Cosmologi- 
cal N-body simulations usually start from initial conditions 
set at some early epoch, then follow the later evolution of 
the dark matter (and, if included, other matter components) 
into the nonlinear regime where it can be compared with the 
large-scale structure in galaxy surveys. In these simulations, 
the dominant mass component, the cold dark matter, is as- 
sumed to be made of elementary particles that interact only 
gravitationally, so the collisionless dark matter fluid can be 
represented by a set of discrete point particles. 

Gravitational, long-range interactions are peculiar be- 
cause: (i) they are inherently negative specific heat sys- 
tems (he;, taking energ y away from it heats it up; see, 
e.g., Lyndcn-Bcll & Lvnden-Bell 1977; Padmanabhan 1990; 
Ihvnden-Belll I2001TI . or, equivalently, internal tem perature 
gradi ents are enhanced instead of being erased fel-Zanti 
1998); (ii) the potential energy is a superextensive quan- 
tity, which leads to an intrinsic instability in the case 
of virialized systems when their mass is increased (e.g., 
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iHeggie fc H ut 2003). The difficulty of characterizing the dy- 
namical evolution and final equilibrium state of collisionless 
systems is deepened by the fact that the classical statisti- 
cal mechanics and thermodynamics have been little devel- 
oped for gravitational problems due to complications posed 
by the long -range, unshielded nature of the gravitational 
force (e.g., |Padmanabharjri990l . 1200a) . Here we focus on 
also peculiar, but largely overlooked property of gravita- 
tional N-body systems, the Energy Ranking Preservation 
phenomenon (her eon, ERP). The ERP was first noted by 
Quinn and Zurek JQuinn fc Z urck 1988); it reemerged later 
in the literature, as a subsidiary result, in different contexts 
such as the evolution of sim ulated galaxy mergers and col- 
lapses JKandrup et al.lll993l) . the numerical study of "vio- 
lent relaxation" processes (e.g., iFunato et al. 1992a), the 
applicability of the sec ondary infall model to collapsing cos- 
mological proto-halos JZaroubi et al.||l 996J). or the study of 
properties of dark matter halos JH enrikse n fc Widrowlll999l : 
iDantas et ai]|2003t iDantas fc Ramosll200d) . ~ 

A gravitational N-body system displays ERP when the 
relative ordination or ranking of sufficiently large groups 
of particles with respect to their mean energy is preserved 
along their possibly intricate trajectories in the phase space. 
In the present context, ERP means that there is a "coarse- 
grained" sense in which the ranking relative to the mean 
one-particle, mechanical energy of given collections of ele- 
mentary particles is strictly not violated during the gravi- 
tation al evolution of the sy stem. As remarked by Kandrup 
et al. (Ka ndrup et alj[l993l) . "mesoscopic constraints" seem 
to be operative at the level of collections of particles when 
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a large N-body system evolves gravitationally towards equi- 
librium. 

The present work is part of a project first motivated 
from a study of the properties of elliptical g alaxies in the 
context of the so-called "fundamental plane" JDantas et alJ 
2003). We have subsequently analysed the ERP in individual 
halos in a cosmological simulation and have found some pre- 
liminary results on the behaviour of the c oarse-grained level 
neces sary for the stablishment of the ERP ( Dantas fc Ramosl 
2006). In the present paper, we aim to address the ERP ef- 
fect within a large-scale N-body cosmological simulation in 
a more general context (as opposed to the restricted "indi- 
vidual halos" approach of previous studies) . We investigate 
the problem in a more global sense (all particles of the sim- 
ulation, inside and outside halos will be considered) and at 
different dynamical regimes and scales, namely, at the linear 
regime, where structures are not yet virialized, and at the 
small-scale nonlinear regime, where particles probably are 
part of virialized or quasi- virialized structures. 

In fact, the non-Hubble component of a galaxy ve- 
locity through the universe (its peculiar velocity), due 
to the acceleration caused by clumps in the matter dis- 
tribution, is an important qu a ntity i n cosmology (e.g., 
Ijuszkiewicz. Sprin gel fe Purred I l999|: J Fukushige fc Sutol 
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Feldman. H. et aZ.ll2003l) . Here we will 
be interested in studying the Pairwise Velocity Distribution 
(PVD) from the point of view of the ERP phenomenon. In 
other words, we analyse a cosmological N-body simulation 
in the linear and nonlinear regimes as a function of different 
energy partitionings, under the perspective of characterizing 
the behaviour of the underlying cosmic flow. 

The present paper is organized as follows: in Section 
2, we present the methodology and illustrate the ERP. In 
Section 3, we analyse the behaviour of the PVD. In Section 
4, we offer a bird's eye view in the literature to contextualize 
the ERP phenomenon within our current understanding of 
the gravitational collapse and virialization mechanisms. In 
Section 5, we present the summary and final conclusions. 



2 METHODOLOGY AND THE ERP 
PHENOMENON 

In the present work, we used in our analysis a A- 
CDM N-body simulation of the VIRGO Consortium 
(the A-CDM model is currently the most observation- 
ally favored cosmological scenario, e.g. ISpergel et alJl2003t 
ISpergel. D. N. et al-lEoOrJ) . 

The simulation box has 239.5 Mpc/h (h = 0.7) of size, 
with 256 3 particles, each with mass of 6.86 x 10 10 Mq/Ij, 
therefore representing approximately galaxy-sized objects. 
We have selected a random fraction of the particles (1 out of 
100) in order that our analysis be computationally feasible. 
We have also restricted our analysis to the initial (z — 10) 
and final (z — 0) outputs of this cosmological run. Parti- 
cles of a given initial model are sorted according to their 
mechanical energies. The energy space is then partitioned 
into a few (e.g. 10) cells or bins of equal number of particles 
and for each of these bins, the mean energy is calculated. 
The mean energies of these same collections of particles are 
then recalculated for the final model and compared with 
their initial values. As we will show below, generally these 
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Figure 1. Top panel: A scatter plot of the initial (rcdshift z = 10) 
versus final (z = 0) energies of the particles in the cosmological 
simulation. Bottom panel: The ERP effect: the ordering of the 
mean energies of collections of particles are mesoscopically main- 
tained for various levels of energy partitionings (Nu na = 10, 20, 
and 50; see text for details). Only the initial and final energy bins 
are shown here, connected by line segments. 



mean energy cells do not cross each other along the evolu- 
tion (ranking is preserved), although particles individually 
gain or lose energy as the system virializes. 

In Fig. we present two panels illustrating the global 
characteristics of the energy space in the Virgo's A-CDM 
simulation. In the top panel, we show a scatter plot of 
the initial (z = 10) versus final (z = 0) energies of the 
particles in the simulation. The mechanical comoving one- 
particle energy of a particle i was calculated classically from 
Ei = l/2mvf — 1/2^ . ,. Gm 2 /|xi — Xj|, with comoving po- 
sition Xi and peculiar velocity Vi. Units used were Mpc for 
length, Gyr for time and Mq for mass. It represents the en- 
ergy with respect to the comoving reference frame (that is, 
the frame which moves with the cosmological expansion of 
the simulation box) . The energy associated to the dynamics 
of expansion of the cosmological box does not enter into the 
above computations. 
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The straight ~ 45° line pattern exhibited in the top 
panel of Fig. highlights the strong memory preservation 
of the initial condition state along the evolution of the sys- 
tem in the energy space. We remark, however, that there is a 
significant spread around the correlation, which indicates in- 
dividual gain/loss of energy ('mixing') during the evolution. 
At the bottom panel of Fig. we present the ERP effect 
for different choices of energy partitioning (Ntins = 10, 20, 
and 50; we have fixed this choice of partitionings through- 
out this work). The energy cells do not cross, meaning that 
the most (least) energetic particles at the initial state are 
still the most (least) energetic particles at the final state. 
Note that the energy ordination is observed despite the 
fact that particles individually have mixed in energy space. 
This unexpected orderly evolution contrasts with the highly 
complex structure of the actual phase space of the system, 
which displays sensitivity to changes in their initial con- 
dition s which are characterist i c of chaotic dynamical sys- 
tems dKandrup fc Smithlll99ll Il992t iGoodman et al.lll993t 



iMerritt fc Vallurll996tlEl-Zanlll99StlKandrup et al.l2003l) 
Hence, clearly, memory is preserved only in a coarse-grained, 
mesoscopic sense (top panel of Fig. . 

Notice also (c.f. top panel of Fig. that some of the 
most bound particles tend to become at the end state even 
more gravitationally bound, as evidenced by a tenuous 'tail' 
at the bottom /left side of the scatter plot (see also Fig. 13 
of | Quinn fc Zurek1ll98SL and Fig. 9 of iHenriksen fc Widrowl 
1999). The final mean coarse-grained energies are dis- 
tributed in a smaller range than their initial counterparts. 
This would not be expected for energies evaluated at a fixed 
background. But in the present case, the energies are calcu- 
lat ed with respect to the co moving frame (for further details, 
see IPantas fc Ram^Einfil . 



3 THE PAIRWISE VELOCITY DISTRIBUTION 

We will proceed our analysis considering the pairwise veloc- 
ity distribution (PVD) from an 'energy space' perspective. 
Here we consider a pair of galaxies lying at a comoving sep- 
aration vector R. We define the pairwise peculiar velocity 



follWs-Cvj-vjJ-R, 



(1) 



where vi and Vj are the proper peculiar velocities of a 
given pair of galaxies in the system. The PVD is the dis- 
tribution of all particles 5v\\ (R) for a given scale R probed. 
The importance of this quantity cannot be underestimated. 
Notice that, separated at cosmological distances, a given 
galaxy pair will be approximately following the uniform 
Hubble flow. In other words, it is expected that the ratio 
5vu(R)/(HR), where H is the Hubble parameter, tends to 
zero in this case. On the other hand, objects at much smaller 
scales evolving under nonlinear regimes will have their dy- 
namics supersede the Hubble flow 1 . 



1 Davis & Peebles JDavis fc Peebles 1977) have argued that for 
an isolated system in virial equilibrium, and therefore highly non- 
linear, the physical (proper) separation of particle pairs within 
the system is expected to remain constant on average (the so- 
called "stable clustering" hypothesis). In the universe, however, 



Hence, the properties of the PVD, evaluated at different 
comoving scales, are intimately related to linear/nonlinear 
regimes. Our major concern here is to describe the behaviour 
of the PVD considering the fact that there is an ordering 
preservation of the mesoscopic particle energy cells during 
the cosmological evolution. How does the nonlinear evolu- 
tion of the pairwise motions behave under such an energy 
constraint? 

Fig. is an illustration {Nuna fixed to 10) of the typi- 
cal results of Skewness and Kurtosis 2 of the PVDs for every 
energy cell separately. All PVDs at z — show significant 
departures from gaussianity. Although this is a behaviour 
well known in unconstrained PVDs, and modelled analyti- 
cally in previous works (e.g., iDiaferio fc Gellerlll996UShetbl 
ll99fAlSeto fc Yokovamall9 98) . it is not obvious why this fea- 
ture should remain in a energy-constrained PVD. In fluid 
turbulence, for instance, energ y-constrained PVDs recover 
gaussianity JNaert et alJll998T) . indicating that the power- 
law tails apparent in the velocity histograms stem directly 
from the spatially intermittent nature of the under lying en- 
ergy cascade that drives the system JFrisch 1995). This is 
clearly not the case here. 

More unexpectedly in Fig. [5] is the evident correlation 
between early and late epoch statistics, for certain scales 
(R=l and 10 Mpc). Note that scales below 10 Mpc are pre- 
cisely those that probe the nonlinear regime, where peculiar 
velocities dominate the Hubble flow and particles are al- 
ready part of virialized or quasi-virialized clusters. On the 
other hand, on larger scales (R ^ 100 Mpc), where there 
is little discernible structure in cosmological box and the 
mass distribution appears homogeneous and isotropic, there 
is no evidence of memory on the skewness and kurtosis of 
the energy-constrained PVDs. 



In Fig. El we plot the dispersion ((8vu(R) 



2 »of 



the PVD as a function of the mean energy cell for the initial 
(z — 10, to the left, circles) and final (z — 0, to the right, tri- 
angles) states. At this point let us review what has actually 
been measured and its meaning: (i) we measure the radial- 
component, peculiar, velocity difference of each pair of par- 
ticles, for each mesoscopic energy bin, to obtain the energy- 
constrained PVD (Eq. [Q. This results in a distribution of 
velocity differences that is well-known to deviate from gaus- 
sianity in the energy-unconstrained cas e (Diafcri o fc Gelleil 
1996; Shcth 1996; S eto fc Yokovamall998l ): and (ii) we mea- 
sure the dispersion of the resulting energy-constrained PVD. 
Notice then that what is being considered it is not the veloc- 
ity dispersion of the particles of a given energy cell. It is the 
dispersion of the PVD for each energy cell. In other words, 
the former and the later are different physical measures. 

Overall, we observe that the more (negative) energetic 
the cell, the larger the corresponding dispersion. Moreover, 
the relation between dispersion and energy follows a power- 
law in which the characteristic exponent increases with the 



no system is exactly isolated, and continuous merging and ac- 
cretion are expected to occur in a hierarchical fashion. Despite 
this complication, the stable clustering hypothesis has been used 
as a model to the power spe ctrum of density fluctu ations encom- 
passing nonlinear structures (Hamilton et al. 1991), and has been 
extensively investigated for halos resulting from cosmological sim- 
ulations (e.g., Jine 2001). 
2 The Kurtosis, as usually defined, was decreased by 3. 
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Figure 2. Skewncss (top figure) and kurtosis (bottom figure) of 
the PVDs obtained separately for every energy cell as a function 
of the mean energy of the cell for z = (solid lines) and z = 10 
(dashed lines). iV(, ins fixed to 10 in all these plots. Dotted lines 
indicate the expectations for a perfectly gaussian distribution. 
The absence of data for z = 10 in the first top panels of the 
figures is due to lack of statistics in this scale (R = 0.1). 



probed scale (c.f. Fig. |3] and corresponding legend). Notice, 
however, that probing smaller scales does not necessarily 
mean that one is also probing more (negative) energetic cells 
of central halo regions, which possess higher velocity disper- 
sions, since, as already mentioned, the velocity dispersion 
and the dispersion of the PVD are different physical mea- 
sures. 

In the case of less (negative) energetic cells (fix the tri- 
angle dots mostly at the bottom of each fitted line in Fig.^J, 
for redshift z — 0, moving to smaller scales in fact decreases 
the dispersion of the PVD (the characteristic exponent de- 
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Figure 3. The coarse-grained mean energies as a function of the 
related dispersion of the PVDs (see text for details). This mosaic 
shows the results for every scale R probed by the estimator; top 
panels: results for an energy partition of Na na = 10; intermedi- 
ate panels: Nf, ins = 20; lower panels: Nf, in3 = 50. Circles: results 
for the initial state (z = 10); triangles: results for the final state 
(z = 0). The absence of points in the first upper panel is due 
to lack of statistics in this scale for z = 10. Straight lines are 
fits to the data points, at z = and iV^ns = 50, disconsidering 
the points corresponding to the boundest energy cells. The cor- 
responding slopes (a) for these fits (according to scale R) are: 
(a, R) = (1.130, 0.1), (1.562, 1), (3.215, 10), (4.167, 100). 



creases for smaller scales). Clearly, these less (negative) en- 
ergetic cells, even at small scales, are probing, most of the 
time, not yet virialized motions, still in the linear regime. 
This r esult is in agreement with those of lDiaferio fc Gellerl 
(1996), where it is found that unrelaxed systems tend to 
have more centrally peaked PVD than gaussian distribu- 
tions. Also, some of the particles composing these less (neg- 
ative) energetic cells, at small scales, can also make part of 
the population of the outermost parts of dark matter halos, 
which are not yet virialized. In any case, all energy cells fol- 
low the same simple scaling law, regardless of the linear or 
nonlinear clustering regime. 



4 VIRILIZATION AND THE ERP 
PHENOMENON 

Up to this point we have shown that the ERP is a phe- 
nomenon that occurs embodying all particles of a A-CDM 
cosmological simulation, regardless of the virialization con- 
dition of the immediate regions to which the particles are 
associated with. Also, it is reasonable to consider that the 
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energy-constrained PVDs scaling laws found in the present 
work are a reflection of a remarkably smooth dynamical evo- 
lution of the large-scale motions resulting from the clustering 
process taking place in the simulation. In this section, we al- 
low for a brief digression and review the relevant literature 
on the virialization of dissipationless gravitational systems 
in order to bring the ERP effect into this context. 



4.1 Dissipationless gravitational systems 

In a first approximation, if the dissipation and evolution 
of the baryonic component of a gravitational system can 
be disregarded, the problem of virialization can be stud- 
ied by the means of numerical experiments that attempt 
to reproduce the purely gravitational evolution of a set of 
N point masses 3 . Well before such experiments could prop- 
erly be perfor med, the theoretical work by Chandrasekhar 
in 1943 (viz., iBinnev fc Tremainelll987l) already indicated 
the existence of different dynamical regimes, namely: if 
(8 In N)/N « 1, the system is collisionless in the sense that 
each particle responds to the average gravitational field of 
the system; otherwise, the system is collisional or marginally 
collisional. 

In the collisional case, local fluctuations or granularities 
of the gravitational potential are important to the overall 
dynamics. The accumulative effect of the energy exchange 
between stars as they deflect each other during an encounter 
results that, at the gravitational two-body relaxation time- 
scale, the "memory" of the initial orbits of the stars have 
been erased. Objects like globular cluster are considered ex- 
amples of such systems. In the collisionless case, however, the 
corresponding relaxation time-scale is longer than the age 
of the universe, hence the system cannot reach a thermal- 
equilibrium state from binary collisions. Therefore, it can 
only relax to a quasi-equilibrium state through some other 
mechanism, which also must operate rapidly enough (see 
detailed discussion and references in the next subsection). 
Elliptical galaxies and dark matter halos are believed to be 
examples of such systems. 



4.2 Brief review on virialization of gravitational 

systems 

The first attempts to provide some statisti- 
cal/thermodynamical description of gravitational stellar 
dynamics and collisionless relaxation started in the 1960s 
with the pioneering work of Antonov (see, e.g., Lvndcn-Bcll 
2001 and refere nces the r ein), Qgorod nikov JOgorodnikovl 
19651). Saslaw JSaslawl Il968t Il969l) . and Lynden-Bell 
IILvnden-Bel]|ll967h . Lynden-Bell attempted to answer the 
problem of how a collisionless system would be able rapidly 
to relax to a quasi-equilibrium state, since the equipartition 
of energy cannot proceed through two-body encounters. 
The resulting theory ("violent relaxation") stated that 
if rapid fluctuations in the mean gravitational potential 
energy field ((</>}) of a system far from equilibrium persisted 
long enough before damping, these oscillations would 
promote the mechanism for changing the energy per unit 



mass (-E*) of each star in the system (dE+fdt = d((j))dt). 
Notice that E+ can change by a large amount under this 
hypothesis, that is, |A_E*| ~ |-E*|, so that the effect would 
be as significant as in the two-body relaxation case. 

Soon after Lynden-BelPs publication, several one- 
dimensional N-body numerical experiments were performed 
to test the outcome of gravitational collapse against the pre- 
dictions of "violent relaxation", but the r esults were mostly 
inconclusive (e.g.. Irlohl fc Febd Il967t iHohl fc Campbel 
19681: iHenor] 119681: iGoldstein et al]ll96ft " 



1969; 



Henonl I1968J: IGoldstein et alJll969i ICuperman et al 



Lecar fc Cohenlll97li see also more recent develop- 



3 Dissipational effects t end to contribut e to the formation of deep 
central potentials (e.g., ICarlberdll986l) . 



ment s along these lines in [Tanckusa 1987J; iMineau et alJ 
1990). These simulations generally started with a given ini- 
tial volume of phase-space uniformly occupied by all parti- 
cles of the system, and the resulting end state object usually 
presented a degenerated core and a distinct h alo population 
containing most of the mass (see also lBouvier fc J anin 1970; 
lGottHl973Tl . which could only be partially fit to Lynden- 
Bell's distribution. These results indicated that "violent re- 
laxation" would probably be an incomplete process, in the 
sense that the fluctuations of the gravitational potential 
would damp too rapidly and not all phase-space cells dynam- 
ically accessible to the system would end up occupied. Con- 
sequently, the final equilibrium product generally resulted 
in a partially relaxed object with a stationary potential (ex- 
cept for small statistical fluctuations) that would not neces- 
sarily represent the configuration that maximized the total 
entropy of the system. 

Notice that the general idea behind most of the ap- 
proaches to derive a theory of collisionless relaxation is to 
assume that the final equilibrium state will be the most 
probable state of the system. Such a state could be pre- 
dicted from first principles by maximizing the entropy of 
the system under the appropriate constraints (usually, the 
total mass and energy), taking into account the fact that 
the phase-space density must evolve under the collisionless 
Boltzmann equation. The standard result turns out to be, 
however, that the entropy is extremized if and only if the 
distribution function is of the isothermal sphere. But the 
energy and mass of an isothermal sphere diverge for large 
radii, so this leads to a contradiction with the initial as- 
sumptions of the problem. The general conclusion is that 
no DF compatible with fixed mass and energy maximize 
the entropy: arbitrarily large entropy configurations can be 
obtai ned by an arrangement of the particles of the system 
(e.g., Bi nnev fc Tremainelll987l) 4 . For instance, the entropy 
can increase without bound by confining a small fraction of 
the mass in a compact core and arranging the remainder in 
a diffuse halo. Although the existence of such a core-halo 
configuration is a natural consequence of the relaxation of 



4 This fact lead to several investigations on the properties of trun- 
cated isothermal spheres, i.e., constrained to a finite spherical shell 
and truncated at a given energy. Although such a procedure seems 
too arbitrary, it can be formally justified observing that the ex- 
ternal parts of a gravitational N-body system behave just like the 
boundary conditions artificially introduced by a fixed shell. Such 
an approach has been proved useful to the physical modelling of 
the problem, since the type of shell introduced (and consequently 
the type of physical quantity that is exchanged with the external 
reservatory) descr ibes the system as a member of a specific kind 
of ensemble (e.g., Padmanabhan 1990). 
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collisional systems, like globular clusters, it does not seem to 
be a physically accessible one to collisionless systems, since 
elliptical galaxies and dark matter halos are not observed to 
have such structures. 

Lynden-Bell was aware that in real systems "violent re- 
laxation" would probably not proceed to completion (see 
section 6 of his 1967 paper). Hence the theory would 
only represent a qualitative idea that needed to be com- 
plemented by other mechanisms and/or superseded by a 
more robust theory. This line of thought was strengthened 
more and more from different fronts, in particular from 
theoretical and observation al studies of t he scaling laws 
in el l iptical galaxies (e.g.. iBinnevl 119 82: Tremameet^dJ 
19861: IWhite fc Naravanl 1987t iRichstone fc Tremainel 198. 



which an explanation in terms of the dynamical evolution 
and expec ted final equilibrium state is still being searched 
for (e.g., iDehnen fc McLaughlir]|2005tlAustinet alJl200^ : 
iBarnes et alJ l2005t iLu et alJ 120051) . Notice that, as pre- 
viously mentioned, simulations of isolated collapses (and 
mergers) tend to produce different end states depending on 
the initial conditions. So the apparent universality of dark 
matter halo profiles in cosmological simulations could be in- 
terpreted as the fact that the underlying relaxation process 
in such scenarios produce systems in a state close to the 
most probable one, as derived according to some more fun- 
damental relaxation principle, yet to be uncovered. 



Trenti et al.l 



Hiorth fc MadserJ|l99lTll99a. Il995l : iTrenti fc Bertinlbooa 4.3 How does the ERP phenomenon fit in? 



an d from mo re el aborate 3D 

body simulations (e.g.. [van Albadal|l9 82: McGlvnn 



ida 

m 



Mav fc van Albadal fl984: IVillumsenf Il984l: lAguilaJ 



Palmer et al]ll99bi iFunato et al.lll992rj) . These simulations 



showed that not all initial conditions would lead to the same 
final end-state; in special, only clumpy, "cold" initial states 
would preferentially lead to objects more closely resembling 
the profiles of elliptical galaxies. 

Further extensions and / or elaborations over 
Lynden-Bell's idea (e.g., _|Severne^^IjUwe3 W8C ; 
Gurzadian fc Sawidvl Il986l : ISpergel fc Hernauistl 11992 : 
Kull et alJ 1997), as well as other developments towar ds a 
more consistent theory of g r avitational relaxation (e.g ., |Shvj 
1 19781: IWiechen et alJll98a IZiegler fc Wiechenl 11983 . | 199C 



[ 19781: 

IChavc 



IChavanis et alJll99tt iNakamural 120001) . were subsequently 
carried out. The main difference between several of these 
approaches lies in the adopted definition of entropy and 
the constraints under which it is maximized. Even though 
there is now a large body of work towards a theory of 
the gravitational relaxation phenomenon, it remains a 
challenging endeavour, and in fact a deeper evaluation of 
some of these theories lead to severe inconsistenci es (viz., 
lArad fc Lvnden-Belllliool lArad fc JohanssorJl2005l) 5 . 

Parallel to the developments is stellar dynamics, N- 
body cosmological simulations of hierarchically clustering 
model universes started to address, at increasingly higher 
spatial and temporal resolutions, the issue of the dynami- 
cal and structural properties of virialized dark matter ha- 
los. For instance, the importance of mergers of smaller sub- 
systems - as opposed to simple collapses - to the forma- 
tion of elliptical galaxies and dark ma tter halos started 
to become clear i mportant issues (e.g., iFrenk et al.lll985i 
IZhang et all2002r) . The mai n result however is the exist ence 
of a universal halo profile JNavarro et al.l ll997. 200^), to 



5 Notice at this point that another physical mechanism that is 
known to contribute to the collisionless relaxation of gravitational 
systems is "phase mixing" , in which the phase-space microscopic 
cells are continually stretched into thinner and thinner filaments 
as relaxation proceeds, until the coarse-grained phase-space den- 
sity no longer changes, and the system is s aid to have reached 
a coa rse-grained steady-state configuration feinnev fc Tremainel 
1987). However, "phase mixing" and "violent relaxation" (what- 
ever the given potentials admit or not the coexistence of regular 
and chaotic orbits) are in fact special cases of a more general relax- 
ation processofcollisionless gravitational systems iZiegler et alJ 
ll994HKandrurJll998ft . 



Clearly, "violent relaxation" is routinely observed in N-body 
experiments, but it rapidly damps out at about the free-fall 
timescale of the system. Gravitational potential fluctuations 
can be violent enough to mix the one-particle phase-space 
structure, but not so violent as to remove all "memory" of 
the initial coarse-grained, mean particle energy ordering. 

Dissipationless relaxation must therefore be described 
as a process in which collections of particles closely spaced 
in energy space can efficiently mix in phase space, but un- 
derlying such a mixing, there is a restriction to the exchange 
of energy, which must proceed within some upper and lower 
bounds. Such bounds can only change linearly and orderly 
with respect to all other energy cells of the system, or at 
least must be restored subsequently to the major poten- 
tial fluctuations. This suggests that a line of approach to 
the EPR effect that might prove fruitful is a formulation 
based on a maximum entropy principle in which a further 
constraint, in addition to the total mass and energy, must 
be imposed, namely, a "mesoscopic constraint", operative 
only at the level of a coarse-grained distribution function 
jKandrup et al.lll993l) . 

Some clues on how this theory must be accomplished, 
at our present stage of understanding, are outlined below: 

• The development of a physically well-motivated, in- 
complete treatment of the "violent relaxation" process will 
be an ind ispensab le s tep to the formulation (see, e.g., 

iHiorth fc Madsenlll99lh . 

• The underlying "mesoscopic" constraint must hold even 
under large fluctuatio ns of the gravitational potential field 
jKandrup et al.lll993l) . 

• Sonre heuristic inputs must be taken into account, 
which we collect from present and previous work, namely: 

(i) evidences that the ERP depends on halo mass, in the 
sense that more massiv e halos show more rank preservation 
than less massive ones JDantas fc Ram os 2006). 

(ii) some reasonable validity of Arnold's theorem when 
applied to individual halos, showing that the "mesoscopic" 
constraint acts at the level where collections of particles 
behave dynamically as an individual particle with a char- 
acteristic frequency (or al ternatively, energy) in the mean 
potential field of the halo JDantas fc Ra mos 2006) . 

(iii) the fact that cosmic flows evolve in an orderly 
sense, when tracked from their coarse-grained energy cells, 
even when nonlinearities are already developed (Figurc|2]of 
present paper). 

(iv) other "memory" -like effects, like those evident in 
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higher-order statistics of the PVDs, for scales R ^ 10 Mpc 
(Figure [5] of present paper). 

(v) other scaling correlations associated with the ERP 
phenomenon, namely, that simulations of galactic mergers 
produce "fundamental- plane" -like objects and some energy 
rank violation; whereas sim ple collapses are h omologous sys- 
tems and show clear ERP foantas et alJl2003l) . 



5 CONCLUSION 

Summarizing, in this paper we have analyzed the behaviour 
of a cosmological N-body simulation in the linear and non- 
linear regimes as a function of different energy partitionings. 
We have found evidences that ERP effect, already observed 
in different, more restricted contexts, is also present in a 
large-scale, realistic cosmological N-body simulation from 
the VIRGO Consortium. We have observed that ERP-like 
"memory" effects are also evident in higher-order statistics 
of the PVDs, for scales R ^ 10 Mpc. Finally, we have noticed 
that there is a general scaling law relating the energy cells 
to the energy-constrained PVD velocity dispersion that ap- 
pears to hold regardless the dynamical regime being probed. 
The expansion of the cosmological background brings 
more complexity to the virialization problem when com- 
pared with the "simpler" studies of statistical mec hanics 
of iso lated gravitating systems (IPadmanabhan fc Engineer! 
1998). There are, however, common effects that persist in 
both scenarios, like the ERP phenomenon here investigated. 
Although several insights can be obtained from the bulk of 
current knowledge on the nonlinear g ravitational evol ution, 
a complete theory is yet to be derived JPadmanabhanl2 005 ) . 
For instance, when viewed from the energy space perspec- 
tive, the cosmic motions and clusterings, at several stages 
of development (from the linear throughout the nonlinear 
phase) , do display a certain degree of simplicity and an un- 
expected level of order. The ERP phenomenon is certainly 
an important ingredient to the outstanding open problem of 
the nature of the collisionless relaxation process. 
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